Anderson localization as a parametric instability of the linear kicked oscillator 
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We rigorously analyse the correspondence between the one-dimensional standard Anderson model 
and a related classical system, the 'kicked oscillator' with noisy frequency. We show that the 
Anderson localization corresponds to a parametric instability of the oscillator, with the localization 
length determined by an increment of the exponential growth of the energy. Analytical expression 
for a weak disorder is obtained, which is valid both inside the energy band and at the band edge. 
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I. INTRODUCTION 



Recently it was shown that quantum one-dimensional 
tight-binding models with any diagonal site-potential can 
be formally represented in terms of a two-dimensional 
Hamiltonian map M. On the other hand, this classical 
map is associated with a linear oscillator subjected to a 
linear force given in the form of time-dependent delta- 
kicks. In this picture, both the frequency of the unper- 
turbed oscillator and the period of the kicks are deter- 
mined by the energy of an eigenstate, and the amplitudes 
of the kicks are denned by the site-potential in the origi- 
nal quantum model. It was shown that by exploring the 
dynamics of this classical system, one can obtain global 
characteristics of the eigenstates, such as the localization 
length defined by the Lyapunov exponent of classical tra- 
jectories. 

Inparticular, analytical estimates have been obtained 
in |ij for a specific diagonal site-potential potential 
with short-range correlations (the so-called dimer model, 
see Q). Other applications to the case of general corre- 
lated diagonal || and off-diagonal disorder have re- 
vealed very important peculiarities. One of the most in- 
teresting results has been obtained in |j| where it was 
shown how to construct random potentials with specific 
two-point correlators which result in the emergence of 
the mobility edge in one-dimensional geometry. Based 
on these predictions, very recently experimental realiza- 
tion of this effect has been done in one-mode microwave 
guides f|. 

In this paper, we perform an analytical study of the 
standard Anderson model with diagonal uncorrelated dis- 
order, paying main attention to the problem of the math- 
ematical correspondence between the quantum model 
and its classical representation in the form of a linear 
kicked oscillator. More specifically, we are interested in 
the connection between the Anderson localization and 
the paramtric instability of the corresponding classical 
system. Although the results obtained for the localiza- 
tion length in the weak disorder limit are well known from 
other studies, the method we use here is a new one and 
it may explain the mechanism of the Anderson transition 



in new terms. Moreover, this approach may be very use- 
ful for 2D and 3D cases, for which analytical results for 
global properties of eigenstates are very restricted. 



II. DEFINITION OF THE MODEL 

In this paper we study the relation existing between 
the standard ID Anderson model and a related physi- 
cal system, a linear oscillator with noisy frequency. The 
quantum model is defined by the stationary Schrodinger 
equation (6) 



1p n+ l + tpn-l + e n 1pn = Etpr, 



(1) 



where ipn represents the electron wave-function at the n- 
th lattice site, and the site-energies e„ are independent 
random variables with a common distribution p(e). In 
the standard Anderson model the probability p(e) has 
the form of a box distribution, 



P(e) 



1 



w 



(2) 



whose width W sets the strength of the disorder. In the 
following, however, we will not restrict our considerations 
to the specific form (||) of the probability distribution, 
but simply assume that the variables e n have zero mean 
((O = 0) and a finite variance (e^). 

The kicked oscillator is a harmonic oscillator that un- 
dergoes periodic and instantaneous variations of the mo- 
mentum (the 'kicks'). Such a system is defined by the 
Hamiltonian 



H 



where 
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£(t) = J2 A n 5(t~nT). 



(3) 



(4) 
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The random coefficients A n that appear in the definition 
of the noise (0) represent the intensity of the 'kicks', i.e., 
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they are proportional to the sudden momentum changes 
experienced by the oscillator at times t = nT. In other 
words, the system ([}]) represents a harmonic oscillator 
with a mean frequency uj perturbed by the noise term 
Using the definition (||), one can easily reconduct 
the statistical properties of the noise to the corre- 
sponding properties of the variables A n ; in particular, 
the mean and the variance of £(t) can be expressed as 
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(£(*)) = lim - / £(t) dt 
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and 



(f(t)£(* + a)) = lim - f 1 i(t)i(t + a) dt 

= lim - Y A 2 n = { -^-S(a). 

n=-N/2 

The equivalence of the models (Q) and (||) has been 
discussed in E|. There it was shown how the two- 
dimensional map 



x n+ i = x n cos (uoT) + (p n - A n x n ) sin (wT) 
Pn+i = —x„ sin (wT) + (p n - A n x n ) cos (wT) 



(5) 



can be derived by integrating over a period T the Hamil- 
tonian equations of motion of the kicked oscillator (|^). 
Note that in the map (|^) x n and p n stand for the coordi- 
nate and momentum of the oscillator immediately before 
the n-th kick. Eliminating the momentum from equa- 
tions (|^), one eventually obtains the relation 

x n +i + Xn-i + A n cos (wT) x n = 2x n cos (ojT) , 

which coincides with the Anderson equation (lj) if one 
identifies the site amplitude ip n with the coordinate Xn 
of the oscillator and if the parameters of the models @) 
and fl|) are related to each other by the equalities 



e„ = An sin (ljT) ; E = 2 cos (uT) . 



(6) 



In Ref. the classical map (||) was used as a tool 
to investigate the properties of the model (0); here we 
focus instead on the direct analysis of the Hamiltonian 
model (||). 



III. THE OSCILLATOR WITH NOISY 
FREQUENCY 

The dynamics of the kicked oscillator (||) is determined 
by the Hamiltonian equations of motion: 



LUp 



(7) 



In order to study the behaviour of the kicked oscillator, it 
is convenient to substitute the couple of differential equa- 
tions (0) with the system of stochastic ltd equations, 



dp = —lux dt 
dx = Lop dt 



x^(AWTdW(t) 



(8) 



where W(t) is a Wiener process with (dW(t)) — and 
(dW(t) 2 ) = dt. 

The systems and (||) can be considered equivalent 
inasmuch as the shot noise £(t) is adequately represented 
by a Wiener process W(t). That is the case if the strength 
of the single kicks is weak, i.e., if the condition 



(Al) « 1 



(9) 



is fulfilled. That can be understood by considering that 
the present situation is analogous to the one that occurs 
in the Brownian motion of a heavy particle suspended in 
a fluid of light molecules. The instantaneous impacts of 
the fluid molecules on the massive particle can be success- 
fully described by a continuous Wiener process, provided 
that each single collision does not produce a significant 
displacement of the heavy particle. When the mass of the 
suspended particle is not much bigger than the one of the 
impinging molecules, the nature of the motion changes 
and the effect of the molecular collisions can no longer 
be depicted by a Wiener process. 

A similar analogy can be drawn between the present 
case and the random walk problem. To be exact, let us 
consider a one-dimensional random walk made by some- 
one that takes steps of length I at times nT (with n in- 
tegral). At each step the walker is supposed to go to the 
right or to the left with equal probability. In this model 
the walker's position changes with each step much in the 
way the momentum of the kicked oscillator does under 
the action of a kick: in both cases the relevant physical 
quantity is varied in a sudden and random way at regular 
time intervals. This analogy makes interesting to notice 
that, by going to the limit 

l->0, T^O, 

while holding fixed the ratio 

T 

the discrete time random walk evolves in a Brownian mo- 
tion with diffusion constant D (see, e.g., Ref. 0). In 
other words, a Wiener process can be regarded as a limit 
case of random walk in the limit of very small and fast- 
spaced steps. In a similar way, the 'jump process' £(t) 
can be described by a 'diffusion process' W(t) when the 
condition is satisfied, with the ratio 



k = 



ujT 



(10) 
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playing a role analogous to that of the diffusion constant 
D. 

This conclusion, although substantially correct, must 
be made more precise. The analogy between the kicked 
oscillator (^|) and the free Brownian particle or the ran- 
dom walker, although of much use, cannot be complete 
because the kicked oscillator, unlike the two latter sys- 
tems, is endowed with an autonomous dynamics dictated 
by the elastic force and independent from the noise. As 
a consequence, in assessing the equivalence of the shot 
noise with the continuous process W(t), one must 
also take into account the possibility that the interplay 
of the discrete noise (||) with the unperturbed motion of 
the oscillator might produce different effects than those 
induced by the addition of the continuous process W(t) 
to the deterministic dynamics of the noiseless oscillator. 
More specifically, one can guess that possible 'resonance' 
effects due to the commensurability of the frequency 1/T 
of the kicks with the frequency lj/2-k of the unperturbed 
oscillator might occur. This is actually the case at the 
band centre, when the two frequencies stand in the ratio 
uiT/2ir = 1/4, as we will discuss in Sec. fv|. Apart from 
this exceptional case, however, the dynamical features of 
the oscillators (||) and (||) do not differ, as will appear in 
the following analysis. 

In conclusion, the kicked oscillator (^) and the stochas- 
tic oscillator (@) are equivalent when the individual kicks 
of the first model are weak. On the other hand, since 
the kicked oscillator (Q) is equivalent to the Anderson 
model (|l|), one can conclude that the stochastic oscilla- 
tor (||) provides a correct description of the Anderson 
model for the weak disorder case. This allows one to 
analyse the solutions of Eq. (||) in terms of phase-space 
'trajectories' of the stochastic oscillator (||). In this pic- 
ture, localized states for the Anderson model correspond 
to unbounded trajectories of the oscillator in the phase 
space, while extended states translate into bounded tra- 
jectories. 

In the following we will restrict our considerations 
to the case defined by the weak disorder condition (|^) 
and thus ensure the equivalence between the Anderson 
model ([l]) and the stochastic oscillator (||). 



IV. LYAPUNOV EXPONENT 

Once we have established the correspondence of the 
Anderson model with the stochastic oscillator (|J) , we can 
proceed to redefine essential features of the first model 
in the dynamical language of the second. In particular, 
we are interested in deriving a formula for the Lyapunov 
exponent, which gives the inverse localization length for 
the eigenstates of the equation (Q). Since these eigen- 
states correspond to trajectories of the stochastic oscilla- 
tor (^) , the Lyapunov exponent is naturally redefined as 
the exponential divergence rate of neighbouring trajecto- 
ries, i.e. through the limit 



A = lim lim — / dt — log — 7^—^ . 

T^oo 8^0 T J 5 6 x(t) 



(11) 



which corresponds to the standard expression 

A' 



1 - 

A = lim — log 



V'n+l 



n=l 
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for the Anderson model (|l|). By taking the limit 5^0 
first, the expression (Thy) can be put in the simpler form 



A = (z) = lim — 



dt z(t), 



(12) 



where the Ricatti variable z = x/x has been introduced 
and the symbol (z) stands for the (time) average of z. 

To compute the Lyapunov exponent, as defined by 
Eq. dl2| ) , it is necessary to analyse the dynamics of the 
variable z = ujp/x. The time evolution of this quantity 
is determined by the Ito stochastic equation 



dz 



(tO 2 + Z 2 ) dt - U]y/(A*)/T dW(t) , (13) 



which can be easily derived using (g) and the standard 
rules of the Ito calculus. 

Notice that, while the position and momentum of the 
oscillator (Q) do not evolve independently from each 
other, the dynamics of their ratio z = cup/x is totally 
autonomous from that of any other variable. As a con- 
sequence, one deals with the single differential equa- 
tion (|l^) instead of having to cope with a set of coupled 
equations like (||). Thus, the introduction of the variable 
z, which is suggested by the definition (|l^) of the Lya- 
punov exponent, turns out to be beneficial also for the 
study of the stochastic oscillator (||) . 

As is known, the Ito stochastic differential equa- 



tion (13) is equivalent to the Fokker-Planck equation \ 7\ 



|m) = |^ 2 +* 2 )mm)] 



u 2 (Ap d 2 p 

2T dz 



(14) 



which gives the time evolution of the probability density 
p(z,t) of the stochastic variable z. In other words, the 
evolution of z(t) dictated by Eq. ([13]) is a diffusion pro- 
cess with a deterministic drift coefficient (uj 2 + z 2 ^ and 
a noise- induced diffusion coefficient lo 2 (A%) /T . 
The stationary solution of Eq. ( |l4| ) is 



p(z) = 



dx exp {$(x/lu)} 



exp{-$(z/w)} 



where G\ and C2 are integration constants and the func- 
tion <&(x) is given by the relation 



§(x) 



x 
"3 



(15) 
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which contains the parameter k defined by Eq. (10). 
Since p(z) is a probability distribution, it must be in- 
tegrable and therefore the constant C\ must vanish. The 
residual constant C2 is determined by the normalization 
condition p(z) dz = 1. The resulting distribution is: 



p(z) 



Nlu 2 



dx exp{$ (x/u) - $ (z/lu)} (16) 



with 



N 




dx 



■. exp 



x 
12 



(17) 



Once the steady-state probability distribution (|16|) is 
known, one can use it to compute the average of z that 
defines the Lyapunov exponent (fL2|) 



A = (z) = 

By this way one obtains 
uj f°° _ 
X= 2N dx 6XP 



dz zp{z). 
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(18) 



(19) 



Formula (]19) is the central result of this paper. It gives 
an expression for the inverse localization length in the 
Anderson model that turns out to be valid both inside 
the energy band and at the band edge (although it fails to 
reproduce the anomaly of the Lyapunov exponent at the 
band centre). The extended validity range of the expres- 
sion (19) is a remarkable feature, because the behaviour 



of the localization length at the band edge is known to 
be anomalous ||[l0| and has been previously derived with 
methods well distinct (and more complicated) than the 
ones used to study the localization length inside the band. 
In the next two sections we will show how expression ( |l9| ) 
reproduces the known formulae for the localization length 
inside the energy band as well as in a neighbourhood of 
the band edge. 

Before proceeding along this line, however, it is inter- 
esting to notice that expressions very similar to those of 
Eqs. (|l6|) and ( pj| ) have been obtained for a different but 
related model: that of a particle in a one-dimensional 
random potential (see, e.g., |ll| and references therein). 
The problem is defined by the continuous Schrodinger 
equation 



ip"(x) + [E-U{x)]$(x) = , 



(20) 



where U(x) represents white-noise, i.e., a delta-correlated 
random potential with zero mean: 

(U(x)) = 0; (U(x)U(x')) = 2DS{x - x'). 

This correspondence is not surprising on two grounds: 
first, Eq. ([y) is the discrete counterpart of the continu- 
ous Schrodin ger equation (|2^) , and second, the station- 
ary equation (|20|) is the formal analogue of the dynamical 
equation for the kicked oscillator 

x(t) + (u 2 +u£(t)) x(t) = 0, 

which represents an equivalent form of the system (|7|). 



V. WEAK DISORDER EXPANSION 



The weak disorder case is defined by the condition (|j). 
This condition implies that, except that at the band edge 
(i.e. for uiT — > 0), the parameter ( |Tc| ) must satisfy the 
requirement |fe| <C 1. In this section we analyse there- 
fore the expansion of the Lyapunov exponent ( |l9| ) in the 
limit k — > 0. This corresponds to studying the behaviour 
of the localization length inside the energy band for the 
Anderson model (|l|) with a weak disorder. 

Making use of expressions (|^) and ( |l7| ) , it is easy to 
verify that for k — > the Lyapunov exponent can be 
written in the form 



A = 



t L(3n + 3/2) 

48 



AT x 

71=0 



(21) 



To the lowest order in k this expression reduces to 

{AD 



A = 



8T ' 



(22) 



which represents the basic approximation for the inverse 
localization length in the weak disorder case. 

Taking into account the relations ([|) between the pa- 
rameters of the Anderson model (|l]) and those of the 
stochastic oscillators (0) and (|[), the variance (A^) that 
appears in formula ( |22| ) can be expressed as 



{Al 



l-E 2 /4 



When the distribution for the random site-energies e n is 
the box distribution (Q), one can further write (e^) = 
W 2 /12-. as a consequence, expression (^2[ ) takes the form 



A = 



1 



W 2 



T96(l -E 2 /4) ' 



(23) 



which, for T = 1, coincides with the well-known stan- 
dard formula for the inverse localization length in the 
Anderson model 

The extra-factor \/T stems from definition (|l2] ) of the 
Lyapunov exponent, which implies that A = (x/x) has 
the dimension of a inverse time. In order to have the cor- 
rect physical dimension, therefore, A must be inversely 
proportional to a time parameter which, on the other 
hand, must be a specific feature of the noise (Q), since 
that is the physical origin of the orbit instability. This 



4 



requirement singles out the period T between two kicks 
as the only parameter which can endow A with the proper 
physical dimension; the proportionality A oc 1/T is thus 
fully justified. 

The expression ( p3| ) corresponds to the result derived 
by Thouless using standard perturbation methods p2[ . 
As such, formula (^3|) fails to reproduce the correct be- 
haviour of the Lyapunov exponent at the band centre, 
where the second order perturbation theory of Thouless 
breaks down and an anomaly appears which was first 
explained in Ref. [^3|. This deviation of the inverse lo- 
calization length from the behaviour predicted by for- 
mula ( p3| ) is a resonance phenomenon, which can be con- 
veniently understood by considering the dynamics of the 
kicked oscillator [lQl . In fact, the band centre corre- 
sponds to the case wT — tt/2 and this equality can be 
interpreted as the condition that the frequency 1/T of 
the kicks be exactly four times the frequency uj/2ir of the 
unperturbed oscillator. This generates a resonance effect 
that manifests itself in a small but clear increase of the 
localization length with respect to the value predicted by 
formula (|2^) for E = 0. Once the origin of the anomaly at 
the band centre is explained in these terms, it is not sur- 
prising that the model (||) fails to reproduce this feature, 
because it is obvious that the Wiener noise W{t) cannot 
conveniently mimic the regularly time-spaced character 
of the shot noise (|]). One might then worry that the 
model (||) provides an inadequate description of the An- 
derson model whenever the period of the unperturbed 
oscillator and that of the kicks stand in any rational ra- 
tio. This is not the case, however, because the resonance 
effect at the band centre is the only one that affects the 
localization length at the second order of the perturba- 
tion theory ( |Hj) and is thus of interest for the present 
work. For all the other 'rational' values of the energy 
E = 2cos7ra with a rational, the effect of the resonance 
on the Lyapunov exponent can be seen only by going 
beyond the second order approximation in the weak dis- 
order expansion (see details in Ref. |Io[|). 



VI. THE NEIGHBOURHOOD OF THE BAND 
EDGE 

Besides the k — > limit considered in the previous sec- 
tion, one can also study the behaviour of the Lyapunov 
exponent (|l9|) in the complementary case k — > oo. Phys- 
ically speaking, the limit k — > oo can be interpreted in 
different ways, depending on the reference model. If one 
bears in mind the kicked oscillator (||), then taking the 
limit k — * oo is tantamount to studying the case of a very 
strong noise. More precisely, the condition k 3> 1 implies 
that the kicks play a predominant role in the oscillator 
dynamics with respect to the elastic force. Notice that 
this is not in contrast with the requirement that each sin- 
gle kick be weak, as established by the condition (^|). In 
fact, regardless of how weak the individual kicks are, their 
collective effect can be arbitrarily enhanced by making 



the interval T between two successive kicks sufficiently 
shorter then the fraction (A^)/u of the period of the un- 
perturbed oscillator. 

From the point of view of the Anderson model (^), the 
analysis of the case k ^> 1 is equivalent to the study of the 
localization length in a neighbourhood of the band edge. 
To understand this point, one should consider that in the 
weak disorder defined by the relation (||) , the only 

way to fulfil the condition k = (A 2 n ) / (ojT) — ► oo is to have 
u>T — > 0. Correspondingly, the energy E = 2cos(wT) 
must approach the limit E — > 2 _ , i.e. the edge of the 
band. Using relations (^), one can also express the con- 
dition k 1 in the significant form 

2-£« (e 2 „) 2 / 3 , 

which shows that the investigation of the case k ^> 1 
corresponds to studying the behaviour of the inverse lo- 
calization length for energy values which are close to the 
band edge on a distance scale set by the fluctuations of 
the random site-potential. 

With the physical meaning of the limit k — > oo clear 
in mind, we can proceed to verify that Eq. ( |l9| ) repro- 
duces the correct behaviour of the Lyapunov exponent 
in a neighbourhood of the band edge. For this, it suffices 
to notice that the substitution 



luT sin (wT) 



k' 



{wry 



(24) 



transforms formula ( |l9| ) in the expression originally ob- 
tained by Derrida and Gardner for the Lyapunov ex- 
ponent at the band edge ||. This implies that Der- 
rida and Gardner's expression coincides with our own 
for ljT — > 0, since in this limit the difference between 
parameters k and kl vanishes. The limit loT — > 0, on 
the other hand, identifies the band edge case: this proves 
that formula ( |l9| ) is correct not only inside the energy 
band (except that for E = 0), but also for E — > 2. The 
extended validity range of expression ([l9]) is a relevant 
and novel feature; indeed, to the best of our knowledge, 
no other formula encompassing the whole energy band 
has been previously found for the Lyapunov exponent in 
the Anderson model. 

To conclude our discussion of the k — > oo limit, we ob- 
serve that in this case it may be appropriate to expand 
the integrals that appear in expressions (]l^) and (|l7|) in 
series of the inverse powers of k. One thus obtains 



A = 



^ (-2^6)" ^ /2n + 3\ . _ 1)2n /3 



4k' 2 



T ~ (-2^6)" (2n + l\ h _ n2n/ 



E 

n=0 



-r 



which is the counterpart of the expansion ( f2l| ) of the pre- 
ceding section. To the lowest order in fc _1 this expression 
reduces to 
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1 s/6 / ujT 



T 2 T(l/6) Vsin(wT) 



2/3 



2V1/3 



so that for _E = 2 (i.e., for wT = 0), the Lyapunov expo- 
nent turns out to be 



A 



1 ^6 

T 2 r(i/6)' 



2x1/3 



in perfect agreement with the result originally found in 
(see also pi). 



VII. DYNAMICS OF ACTION-ANGLE 
VARIABLES 

To gain further insight in the dynamics of the stochas- 
tic oscillator (||), it is useful to consider its description 
in terms of action-angle variables. To this end, we intro- 
duce the polar coordinates (r, 9) through the standard 
relations 



[0 : 2ir]. Things are different when ui — > 0, i.e., at the 
band edge, for in this case both the drift and the dif- 
fusion coefficient tend to vanish in a neighbourhood of 
9 = and 9 = tt. Therefore, for uj — > 0, the projection 
of the representative point onto the unit circle tends to 
stick at the critical points 9 = and 9 = tt; as a conse- 
quence, one can expect the invariant measure p{6) to be 
sharply peaked around these two values of 9 and almost 
zero everywhere else. This peculiar behaviour for lu — > 
gives an intuitive feeling for the origin of the anomalous 
scaling of the Lyapunov exponent with the noise strength 
at the band edge. 

The stationary solutions of (|2^) are periodic in 9 with 
period tt. In the interval [0 : tt] the solution p(9) which is 
bounded as 9 — > tt and is normalized with the condition 
Jo ^ p{@) d0 = 1 can be expressed in the integral form 



P(0) 



1 



2N sin 2 9 



dx 
sin 2 x 



exp {$ (cot x) - $ (cot 9)} 



(29) 



r = \/ x 1 + p 2 
9 = arctan (x/p) 



(25) 



Starting from this definition and Eqs. (||), it is straight- 
forward to obtain the following pair of stochastic Ito 
equations 



dr = ^-r sin 4 1 



2T 



d9 = 



T 



dt- 



cos f sin 



T 



•sin 9 cos 9 dW(t) 



(26) 



dt - 



(A^ 
T 



where $ is the function ( |l5| ) and N is the normalization 
constant (0); outside the interval [0 : tt], the proba- 
bility p(9) is defined through the periodicity condition 
p{9 + kir) = p(9) with k integer. 

Notice that at the points 9 = and 9 = tt the invari- 
ant measure p(9) takes the value p(0) = p(ir) = k/4N. 
Expanding expression ( |l7| ) in power series of k and 
it is possible to show that in the opposite limits k — > 
and k — > oo one has 



dW(t) (27) 



p(0) = p(tt) 



k 



A simple examination of Eqs. ([26|) and (J27J) reveals that, 
while the action variable is coupled to the angular one, 
the latter evolves in a independent way. Therefore, one 
can determine the probability distribution for the vari- 
able 9 by solving the Fokkcr-Planck equation associated 
to the Ito equation @ [§§], i.e., 



4iV(fc) 



3 1.1/3 , f) (i /U2/S\ 

6 1 /6(8 7r )V2r(l/6) K +^^ 1 /«- ) 



for k 
for fc 





oo 



dp 

at 



(9,t) 



d_ 

89 
{AD d 2 
2T 89 2 



T 



cos v sm 



P (9,t) 



[(smH) p{9,t)} 



(28) 



In geometrical terms, the introduction of the angular 
variable 9 is equivalent to the projection of the point 
(x(t),p(t)), representative of the stochastic oscillator (||) 
in the phase space, onto the unit circle by the rela- 
tion (p5|). A glance at the Fokker-Planck equation ( |2"g| ) 
shows that the ensuing motion of the projected point 9(t) 
is the combination of a drift with a noise-induced diffu- 
sion. In the Fokker-Planck equation (^8|) the diffusion 
term vanishes for 9 = and 9 = tt, but the drift coeffi- 
cient is not zero at these points (unless uj — > 0), so that 
it is reasonable to expect the invariant measure p(9) to 
be finite and non-vanishing at any point of the interval 



These expansions show that, upon increasing k, the prob- 
ability distribution (|2^) does become more and more 
sharply peaked at the points 9 = and 9 = ix\ the anal- 
ysis of the function p(ff) thus confirms the singular be- 
haviour of the invariant measure that was suggested from 
physical considerations on the dynamics of the stochastic 
variable 9(t). 

Using (^9|) one can easily compute the Lyapunov ex- 
ponent as defined by the expression (|l2|). To this end, 
it suffices to notice that in polar coordinates the Ricatti 
variable takes the form z = u> cot 9: one can therefore 
evaluate the Lyapunov exponent as 

A = w(cot 9) = 2u f d9p (9) cot 9 = 
Jo 

' cot 9 exp{-$ (cot#)} 



exp {$ (cot x)} 



N j Q sm' 
w dx 



■ 2 

sm x 



G 



With a chan ge o f variables this expression can be reduced 
to the form (Jig), so that one eventually recovers the re- 
sult (|l^) for the Lyapunov exponent. This is not sur- 
prising, since the distributions ( |l6|) and (|2^) are strictly 
related; indeed one has 



p(9) = - P (z(e)) 



dz 
IB 



where the factor 1/2 derives from the normalization con- 
ditions chosen for the two distributions. From a geo- 
metrical point of view, the passage from the process z(t) 
to 9(t) is equivalent to projecting the real line — oo < 
z(t) < +00 onto the unit circle < 0(t) < ir via the 
relation 9(t) = cot -1 (z/ui). 
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